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ABSTRACT 

Determining radiation location observationally plays a very important role in testing 
the pulsar radiation models. One-photon pair production in the strong magnetic field, 
7 — e^e^, is one of the important physical processes in pulsar radiation mechanisms. 
Photons near pulsar surface with sufficient energy will be absorbed in the magneto- 
sphere and the absorption optical depth for these GeV 7-ray photons is usually large. 
In this paper, we include the aberrational, rotational and general relativistic effects 
and calculate the 7-B optical depth for 7-ray photons. Then we use the derived op- 
tical depth to determine the radiation altitude lower bounds for photons with given 
energies. As a case study, we calculate the lower bounds of radiation altitudes of Crab 
pulsar for photons with energy from 5 GeV to 1 TeV. 
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1 INTRODUCTION 

Pulsars are ideal laboratories to test physics rules at extreme environment, yet how the pulsar radiations are generated in pulsar 
magnetospheres is surely an open question though more observational data are accumulated. Three t ypes of models have been 
prop o sed to explain the pulsed 7-r ay radiation from pulsars. These models are the p olar cap model ( Ruderman fc Sutherland! 
I975I: Daughertv fc Hardine 1994 ) , the annular gap model joiao et ah 2004 . 2007 ). and the outer gap model jcheng et al.l 
1986al lbl: IChiang fc Romanil] 19941 : ICheng et al.|[200ol : I Zhang et aflboOTl ). The distinguishing characteristics of these models 



are the acceleration region for primary particles. The primary acceleration regions locate near the pulsar surface for polar 
cap models and for the annular gap model Q but far away in the magnetosphere for outer. The classical outer gap interacts 
with the annular g ap due to the p article out flow out from the annular gap, which behavior similar to the outer gap with 
boundary currents ( HirotanTlbood ) . The fundamental physics differs in three types of models due to the difference of magnetic 
field intensity. Observational judgments for the radiation altitudes are therefore very critical to d iscriminate these radiation 
models. Early observations indicated a 10 GeV cut-offs in the photon spectrum of 7-ray pulsars ( Thompson 20081 ) . Thanks 
to the improvements of 7-ray observational technologies, recent results both from space-based Fermi Large Area Telescope 
(FERMI) and from ground-based Cherenkov telescop e (MAGIC) confirm such spectra cut-offs and also have observed the 
high energy tails of the spectra cut-off above 10 GeV ( Lopez et al. 20091 : Fermi Collaboration 2009al lbl). These existing high 
energy observations shied lights on the altitudes of the pulsar 7-ray radiation regions. 



The first order relativistic quantum processes allow interchange of momentum between photons and magnetic field (jErber 

19661) . The electron-positron pair production becomes no longer forbidden, such that a 7-ray photon with sufficient energy can 
convert into e^e~ pairs in strong magnetic field. The cross section for the 7-B process increase with the background magnetic 
field and photon energy. Since the magnetic field near pulsar surface is usually strong, the 7-B absorption process is capable 
of developing a spectra cut-off in high energy ends of the pulsar spectra, if 7-ray photons are generated in the vicinity of star 

surface. 

HardingI (|200ll ) and iBaring et al] (| 19971 ). iBarinj (|2004l ) find an analytical 



After summarizing the works of IBaring 



formula to calculate such spectrum cut-off energy as a function of pulsar period, the pulsar surface magnetic field and the 
radiation loc ation alt i tude. One can fit the observed 7-ray spectra to a theoretical model to measure the spectra cut-off energy, 
and then use iBaringI |2004l )'s formula to determine the lower bound for radiation altitude. However in order to measure the 



^ The annular gap model has extended radiation region. 
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spect rum cut-off energy, one need to include both the absorption effects and the initial photon source function ([Baring et al 



19971 ). which is sensitive to the input radiation model. Naturally one could get different values for the cut-off energy for the 



same data set by using different theoretical models. So determining the radiation altitude lower bounds by using spectrum 
cut-off energy is rather model-dependent, becase there is no model independent way of defining the value of cut-off energy. 

In this paper, we present a model independent method to measure the radiation altitude lower bounds. In fact, one can 
answer the question about the lower bounds of radiation altitude through a slightly different approach. Instead of using cut-off 
energy, one can ask what is the lower bound for radiation altitude, such that 7-ray photons with certain energy e will not 
be absor bed and can be observed? H ere, the term 'not absorbed photons' and 'observable photons' means those optical-thin 
photons ( Rvbicki fc Lightman 1986l l. i.e. photons with optical depth r less than a prescribed optical depth threshold rth, 
where rth is used to specify the confidence levels for low bounds of the radiation altitude. It turns out that there exists a 
minimal radius (called last trapping radius, rti), below which the photons are mostly absorbed and convert into e^e~ pairs 
in the presence of intense magnetic fields. The existence of such last trapping radius is a consequence of the r~'^ dependence 
of dipole magnetic field of pulsar, i.e. we can find a small enough r such that the magnetic field is strong enough to initiate 
the 7-B absorption for the 7-ray photons with sufficient energy. When pulsed 7-ray emission from a pulsar is detected at a 
given energy band we can determine the rtr correspondingly. Obviously such rtr is the lower bound for radiation altitude for 
photons with certain energy. Since such lower bounds only involves the optical depth for absorption processes and observing 
energies, the rtr is a model independent altitude lower bound. 

In this paper, we dedicate to investigate the 7-B absorption radiation transfer in the pulsar magnetosphere. The rtr is 
calculated as a function of the magnetic field configuration of pulsar magnetosphere and the photon energy. Analytical formula 
are presented in i) l2.1l The more precise numerical calculations are made in i; 12.21 for different source locations and in 12.31 for 
phase resolved absorption. The application to Crab pulsar with the phase resolved lower bounds for radiation altitudes are 
presented in § 12.41 Conclusions and discussions are given in § [S] 



2 7-B ABSORPTIONS OF HIGH ENERGY PHOTONS IN PULSAR MAGNETOSPHERE 
2.1 7-B absorption — the analytical approach 

In this section, we analytically estimate the optical depth and the radius for last trapping surface rtr. The rtr is mathematically 
defined as the radius, at which the 7-B absorption optical depth is equal to a given threshold optical depth rth. The probability 
is less than e"'^"' for photons to escape from a regi on with altitude lower than rtr. For examples, when rth ~ 1, 10, the escape 
probabilities are 37% and 5 x 10~^ respectively ( Rvbicki fc Lightman 19861 ). Obviously, such probability tells about the 



confidence of rth as an altitude lower bound. 

Here we only consider magnetic field much weaker than 10^^ Gauss, and the photon energy is much larger than e'^e 
static mass, i.e. much larger than 1 MeV. In this way, we can ignore the photon splitting effect ( Baring fc Hardin3 2001 



Mean while the number of pair stat es is large enough to allow us ignore the resonance effects due to the finite number of pair 
state s I Daughertv fc Harding 19831 ). The following formula is appropriate to calculate the 7-B absorption coefficient k ( Erbei 



1966h 



K = 1.55 X 10^efM%-^V3 ) cm-V (1) 

The function A'l/a is the second type modified Bessel function of order 1/3; -Bx[i2] is the magnetic field strength perpendicular 
to the photon propagating direction in unit of 10^^ Gauss; and e[McV] is the energy of photon in unit of MeV. 

The photons of 5 GeV to 1 TeV in the pulsar magnetosphere are radiated by ultra-relativistic charged particles, which 
are approximately beamed along the magnetic field lines. After a photon propagates a length of A, the angle 9i between 
photon direction and local magnetic field becomes di ^ A/rcur. Here rcur is the curvature radius of magnetic field lines. 
For dipole magnetic field, we have rcur ~ t^^^^^ic'^' where r is the altitude of the photon and ric is the light cylinder 
radius. Using geometrical unit such that light velocity is 1, we have ric ~ p/{2n) for the pulsar with period of p. One gets 
B_i_ ~ 6iB — A/rcur-B, where B is the magnetic field intensity at the altitude. The dipole magnetic field has intensity of 
B ~ BoB? /r^ , with Bo as the magnetic field intensity at the pulsar surface, and R as the radius of the pulsar. The optical 
depth r = KdX. Because the absorption coefficient k exponentially dependent on the magnetic field, which is of r 
dependence, the absorption coefficient exponentially decrease as e^^ when A is large. The dominant part of the optical 
depth are thus the integration of A from to rs, where rs is the altitude, at which the photons are generated (see the Appendix 
for the details). Further more, when A ^ rs, the magnetic field B can be approximately regarded as a constant. Thus the 
characteristic absorption length Ac is Ac — rs, and the integration for optical depth is then replaced by r ~ kAc ~ Krs, which 
is later further justified by numerical integration. Thus 
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Figure 1. The radius of last trapping surface rtr for different photon energy e and surface magnetic field. The x-axis is the photon energy, 
while the star surface magnetic field intensity Bo are labeled above each curve. For each Bq there are two groups of curves. Two solid 
curves are the numerical result from directly solving T(rtr) = t^jj with Tth = 1, 10 from the top down. Three dashed curves are analytical 
solutions from Eq. UJwith rth = 1, 10. 



e[MoV] \^ -D0,[12]-K[cjj,]e[McV] J 

where rs_[cm] is the Vs in unit of cm, is the period of a pulsar in unit of second. 

By setting T(rtr) = Tth and solving for rtr, we get radius of last trapping surface, which should be the lower bounds for 
the radiation altitude for photons. The numerical solutions to T(rtr) = Tth is given in Fig. [T] for diflferent Tth- 

We can also get analytical solutions to T(rtr) = Tth using asymptotic method. Note that the Bessel function have 
asymptotic property of K^^^{x) ~ ne''^^ /{2x). The solution to Eq. [2]can be presented by usi ng the Lambert W function, 



where function y — W(x) is defined as the solution to equation ye^ = x i Corless et al. 19961 ). Using the function W and 
taking R — 10® cm, the solution to T(rtr) ~ Tth is 



-1/5^2/5 / 3.46 X 10^°<,^,; 

^,[12]e[MoV]P[s] 1/3 5/3 



(3) 



which have asymptotic approximation given as 

nr.M ^ (1.1 X 10^ - 1.9 X 10^ InrtOBo./^^jei^tvjPH (4) 

This analytical results are plotted together with the numerical results in Fig. [1] The analytical approximation give error of 
about 20% percent. The effect of Tth is very tiny, due to the logarithmic dependence of Tth as we expect from Eq. 2] We 
present the results with Tth = 1, 10. It turns out that the difference between the rtr for Tth = 1 and for Tth = 100 is still less 
than 20%, where the probability for photons escaping region with Tth = 100 is 10^*'^ times smaller than escaping region with 
Tth = 1, or the radiation from region with Tth = 100 is 10'*"^ timeas fainter than radiation from Tth = 1- This is far beyond 
dynamical range of instruments and is exactly why we can define the optical thin photons as observable photons and use rtr 
as the altitude low bound, without going into the details of equi pment response. 

For the case of Crab pulsar, i.e. Bo, [12] = 3.7, p[s] = 0.033 ( Manchester et al. 20051 ) . we see from Fig. [T] that the radius 
of the last trapping surface for 5 GeV, 10 GeV, 100 GeV and 1 TeV 7-ray radiation are 100 km, 150 km, 350 km, 850 km 
respectively. 



2.2 7-B absorption processes with detailed geometry 

We have estimated the radius of last trapping surface rtr and taken it as lower bounds of the radiation altitudes for the pulsed 
pulsar 7-ray radiation. The 7-B process dependents on geometries of both photon ray directions and magnetic fields. In this 
section, we calculate the rtr with detailed geometry. 

The 7-B absorption coefficient k is sensitive to the geometrical parameter 9i, the angle between the photon propagating 
direction and the local magnetic field direction. There are five possible geometrical effects that changes the 6i. These effects are 
listed as follows. (1). Photons, when generated, are not exactly beamed along local magnetic field lines du e to the finite Lorentz 
facto r 7p of particles. This introduces an increase of Oi by SOi, where Sdi ~ 7p ^ < 10"^ given 7p ^ 10^ (jRvbicki fc LightmanI 



198' 



3)- (2). Due to special relativistic effect, the rotational veloci ty of the pulsar bo osts the photons to a new direction, which 



gives SOi ~ J7r ~ r/ric ~ 10~^ for a young pulsar with p ~ 0.1 s jGangadharalboOsI ) . where Q, — 2tv/p is the pulsar's rotational 



angular velocity. The r is choose to be 100 km as order of magnitude estimation from previous section. It should be noted 
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Figure 2. The geometrical configuration for the coordinate systems for calculations. In the laboratory frame O — XYZ, the Z axis is 
aligned with the rotational axis CI of pulsar. The magnetic frame O — X'YZ' is generated by rotating the O — XYZ coordinate around Y 
axis by an inclination angle of a, such that the Z' axis is aligned with the dipole magnetic momentum /i. The polar angle and azimuthal 
angle in laboratory polar coordinate are denote as 0, 3>, while the polar angle and azimuthal angle in the magnetic polar coordinate are 
6, ip. The rg is the radius for the radiation source. 



that we use geometrical unit through out this paper without special mention, i.e. the hght speed c = 1 and the gravitational 
constant G = 1. (3). The magnetic field is co-rotating with the star, thus after the photon propagates for time of t, the 
magnetic field has rotated for an angle of Qt. It gives 59i ~ f2£ ~ 2Tvr/p ~ 10~^ . (4). Due to the gravitational field of pulsar, 
the photons propagate along 'curved' geodesies rather than 'straight' lines. This efi'ect is order of 56i ~ dirajr ~ 10^"^, as 
m ~ 2 km for 1.4 solar mass pulsar. (5). The pulsar drags the background space co-rotating with it due to gravitomagnetism 
effects. The space co-rotating angular velocity is Lense-Thirring angular velocity cjlt ~ 2mQ,R^ /r"^ , where R is the pulsar 
radius. This gives 59\ ~ uJurt ~ 2niQ,Ft^ /r^ ~ 10"'^. We see that the curved spacetime effect and frame-dragging effect s are o f 



higher order geometrical effect compared to other effects, which are coincident with the results of lGonthier fc Harding (| 19941 ). 
Therefore, we only consider the aberration effect and the magnetic field rotation effect to correct the geometrical effects to the 
order of 10~^. This conclusion are valid for short period pulsar (p ^O.ls) and higher altitude absorption (r ^lOOkm). Thus, 
we can just use flat space-time geometry (including photon direction and magnetic field direction) to calculate the geometrical 
parameters (e.g. ft). 

However, the gravitational effects also play the other two important roles in 7-B processes. (1). The gravitational field 
makes the photon red shifted, i.e. the energy e of photons observed by static observer is e = £00(1 + ??), where Eoo is the photon 
energy observed by infinite-distant observer and rj = m/r is the Newtonian gravitational potential. This i ntroduces energy cor 



rection by em/r ~ 10 e. (2). The dipole magnetic field strength is enhanced due to grav itational effect IjGinzburg fc Ozernov 



1964 Underson fc Cohenlll970l : IWasserman fc Shapirolll983l : iMuslimov fc Tsvganlll986l l. which magnetic field intensity cor 



rection of order m/r ~ 10"'^. Both of the two effects are of 10~^ order effects. In short, we can use flat space geometry 
to calculate the geometrical parameters, i.e. 6i, but the magnetic field and photon energy need the first order gravitational 
corrections to include effects of 10"'^ level. 

We set up the coordinates as shown in Fig. [2] In the coordinate O — XYZ, the Z axis is aligned with the rotational 
axis n of pulsar. This coordinate O — XYZ is called the laboratory frame. Another coordinate, the O — X'YZ' , is generated 
by rotating the O — XYZ coordinate around Y axis by inclination angle a, such that the Z' axis is aligned with the dipole 
magnetic momentum ^. We call this coordinate the magnetic frame. The two vectors ft and /i locates in the plane O — XZ, 
which is called — jj, plane. The polar coordinate associated with O — XYZ and O — X'YZ' are called the laboratory polar 
coordinate and the magnetic polar coordinate, respectively. The polar angle and azimuthal angle in the laboratory polar 
coordinate are denoted as 0,$, while the polar angle and azimuthal angle in the magnetic polar coordinate as 9,ip. Here, 
we use bold type font to label the vector or the matrix; while we use subscript x.y,z and x' ,y' ,zi to indicate their component 
value in the laboratory frame and component value in the magnetic coordinate respectively. The coordinate transformation 
between two coordinates is given by a' = T'-' • a'-", where a and a' are any vectors in laboratory frame and magnetic frame 
respectively, transformation is the matrix given as 

cos OL sin OL 

TL' = ( 1 I (5) 

- sin OL cos a 



Let's consider a photon generated at the position 'S' in Fig. (2] The initial position of photon Ts in the magnetic polar 
coordinate is 
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sin 6s cos ifs 
sin 6s sin (ps 
cos 6= 



(6) 



where 6s, ifs are the polar angle and azimuthal angle for source position 'S' in magnetic frame, and rs is the radiation altitude. 
The photon generated at position 'S' propagates along direction hb of the local magnetic field fo r a co-rotating observer. The 
component value of vector ne is easily calculated in magnetic polar coordinates ( Swisdak 20061 ) 



?^B,z' 



V 



3 cos (?g cos ips sin (?g 

^l + 3cos2 
S cos 9s sin sin 0s 

yi + 3cos2 
1 + 3 COS 2^3 



-^1 + 3 cos^ I 

which can be used to calculate the component value of ns in laboratory frame by combining it with Eq. [S] 

The rotational effect boosts the photon to direction n^, k nown as the aberration effect. The relation between components 



of riiy and components ub is give by Lorentz transformation IjMisner et al.iil973i ) 
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where matrix C is defined as 

7«i^ l + (7-l)^ 
(7 - 1 
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(7-1)- 

1 + (7 - : 









iij — and P = + Vy respectively, x is the redshift of photon due to aberration effect. The 



and 7, /? are 7 = — "'-^ 

Wa; , Vy is the components of co-rotation velocity at the place where photon is generated. The velocity can be calculated by 
V = X Ts. In the laboratory frame, we have 

— sin ^s sin (fs 
cos ds sin a + cos a cos Lps sin 6s 


After take such aberration effect into account, the equation for component value for is 









\ =rsni 







3/i + / sin a — rsO, sin 0s sin ips [^/2b-\-3drs^ sin(2^s) sin ips] 

V2nbx 

n a.)] sin 0^ sin 



/ COS a — 3i 



2gbx 

ih tan a 



(7) 



V2bx 



where 



9 = 



1 — rs |^((;os a sin 9s cos yJs -|- sin a cos &s)^ -I- sin^ 9s sin^ y>s] , 

6= ^3cos(26is) + 5, 
d = cos Q sin 6 cos (/9s -I- sin a cos 6's 
/ = 3cos(26»s)-(-l, 
h = cosasin(20s) cosies, 

sin Q! sin 0^ sin tps 



The position of photon x(t) is a function of photon propagating time t. With the n^, we can calculate the position of 
photon by x(t) — Ts + n^t. We need to further calculate the magnetic field at each position x{t) to get the absorption coefficient 
K. To do this, we firstly calculate magnetic field in the magnetic polar coordinate and then use coordinate transformation to 
derive the magnetic field strength at x. 

The magnetic field in the co-rotating magnetic polar coordinate is 



B' 



2r» 



1-1-?? 



8 + 4 cos(26l) \ , / 5 + 3 cos(26l) 
5-|-3cos(26l) 



-na 



(8) 



where T,9,ip are the radius, polar, and azimuthal angle in the co-rotating magnetic polar coordinate for x(t). Bo is the 
effective surface magnetic field determined from energy losing rate, and 77 = m/r is the Newtonian gravitational potential. 
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The magnetic field configuration defined in Eq. |8] satisfy two conditions that the magnetic field intensity is concord with 
results in the Schwarzschild background and the magnetic field direction is concord with the flat background results. Due 
to the rotation of pulsar, the magnetic field is also 'rotating'. Then the relation between magnetic field in laboratory frame 
and the magnetic field in the co- rotation magnetic polar coordinates is B = Tn ■ Ta ■ B', where the transformation Tn is 
coordinates transformation due to the rotation of the pulsar. The matrix form of Tn is given by 

cos Q,t — sin Q,t 
Tn = I sin D,t cos D,t 
1 

Use this transformation, we can also figure out the relation between the co-rotational magnetic polar coordinates {r, 6, (f} and 
laboratory frame {x,y,z}, which is given as 

cos 6 sin a + cos a cos (f sin 6 

sin 9 sin (p \ . (9) 

cos a cos 6 " cos (p sin a sin ( 

Since the photo propagating direction does not change in laboratory frame, the perpendicular magnetic field strength _Bj_ is 
B± — ^/B ■ B (B ■ n^)-^. When the gravitational redshift of photon is taken into account, we have e = £cx)(l + v)- In this 
way, the absorption coefficient k as a function of t can be calculated using the e, B±, and Eq. [1] We then integrate k over t 
to calculate the optical depth. 

We now summarize the steps to calculate the optical depth and the radius of last trapping surface, as follows, 

(i) For a radiation source located at {ra,9s,ips}, we calculate the photon propagating direction Ui, with Eq. [T] 

(ii) We calculate the photon position after time t using x(t) = rs + n,^t; 

(iii) Solve Eq. |9]to calculate the coordinate of photon in magnetic polar coordinate {r,6,if ) from coordinate of x(t); 

(iv) Using the photon position(r, 0,ip ) to calculate the magnetic field strength according to Eq. [8l Then calculate B± using 
Bi_ = ^B-B-(B-n„)2; 

(v) Using r(rs) — K{t\rs)dt to calculate the optical depth for the photon coming from position rs. The k is calculated 
from Eq. [H 

(vi) Solve T(rtr) = nh respected to rtr to determine the lower bounds for radiation altitude given required rth; 

There is no analytical solution to T(r-tr) = 1 respected rtr- We use bi-section method to solve it numerically, while the 
integration is performed using adaptive integration method to refine a preset logarithmic mesh of t to achieve necessary 
numerical precision. The results are given in Fig. [3] We also get the results for rth = 3 and rth = 10, the plots are very similar 
to Fig. [3] due to the logarithmic dependence of rth- 

The results are different for various inclination angles of pulsars, due to the aberration effect and rotational magnetic field 
effect. If these two effects would be ignored, the results would be insensitive to inclination angle a. We also see the aberration 
effects become important for large inclination angles, since the rs shows larger difference between tp^ = 90° (for the leading 
photon rays) and ips = 270° (for the trailing photon rays) cases. 



2.3 Radiation Geometry and Phase- resolved Lower Bounds for Radiaiton Altitudes 



In this s ection we determ i ne the radiation locat i on for different longitud inal phase. The details of radiation geometry can be 
found in Icil et all (|l984l ): lUvne fc Manchester] (|l988l '): iLee et all l|2006l ). We omit the aberration effects here, because it is 
second order effect for calculating the pulse phase. 

Given the pulse profile longitu de A$ (see Fig. [2] for details) and the view angle C,, the half angular beam width for the 
radiation beam 6^^ can be solved bv( Gil et al J 1984: .Lvrie fc Manchester,, 198& ) 



(I) 



M 



V 2 / 



sin a sin C, + sin 



(10) 



With the 9^, the radiation source polar angle 9^ and azimuthal angle tp^ in magnetic polar coordinate can be calculated, 
where 



yjs\n^9^ - lOsin^ 61^ + 9 - sin^ 6I„ 



(11) 



and 
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Figure 3. The lower bounds rtr of altitude as functions of radiation source polar angle 6s. The solid lines, dashed lines and dotted lines 
correspond to ipe = 0°, (pa = 90°, and <ps = 270° respectively. The surface magnetic field strength is taken to be So = 3.7 X 10^^ Gauss 
as the Crab pulsar, and the threshold optical depth is set to = 1. The e is the photon energy used in calculations for each panel, 
while the inclination angle of pulsar a is labeled on the top of each column. The results for Ttjj = 3 and rtji = 10 are almost the same as 
results presented here. 



Using Eg. 1101 [TTl and 1121 we can calculate the angular position for radiation source 9s,tfis from pulse profile phase A$ 
and view angle The we can determined the rtr using the techniques developed in § 12.21 For — 60° case, the result could 
be found in Fig. O 



2.4 Application to Crab pulsar 

As a case study, we apply the method for altitude lower bounds to Crab pulsar and obtain the phase resolved altitude lower 
bounds f or 25 GeV 7-ray emission. T he vie w angle of Crab pulsar is taken to be 60°, as measured from X-ray image of Crab 
nebular JWeisskopf et all [20 00: Wan3 l2003l ) and a ~ 55° — 60° by fitting slot gap model or a ~ 70° from fitting outer gap 



model I Fermi Collaborationll2009bl ). The surface magnetic field of Crab is taken to be Bo ~ 3.7 x 10^^ Gauss. The inclination 



angle for Crab is unknown at this time, thus we calculate the rtr for the four cases, a = 15°, 45°, 75°, 86° 0. Assuming the 
radiation are from single pole, we find that the phase resolved lower bounds are roughly around 250 km for various inclination 
angles, which is concord with our estimation in i\2.1\ 

For pulsar with inclination angle a and the view angle we can use Eg. 1101 [TT] and [12] to calculate the magnetic polar 
angle 6s and magnetic azimuthal angle 0s for radiation source, which contribute to pulse profile with longitude A<I>. The we 
use bi-section method to solve T(rtr) ~ nh respected to rtr to determine the altitude lower bounds, below which the photon 
escape probability is less than e"'^"'. The numerical results are given in Fig. 3] for rth ~ 1 case. 



3 DISCUSSIONS AND CONCLUSIONS 

In this paper, we derive the basic formalism for calculating the effect of 7-B absorption for 7-ray photons in pulsar magneto- 
sphere. We have considered magnetic field line bending, aberration effect, rotational effects, gravitational red-shift of photons 
and gravitational effect on magnetic field intensity. 



2 iRankinI lll990l) got a = 86° 
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Figure 4. The altitude lower bounds for different pulse longitudes. Each column corresponds to results for different inclination angle 
a, as labeled on the top. The x- axes are the puls e longitudinal phase. The first row is the pulse profile of Crab observed in above 10 
MeV band, the data comes from llThompsonll200ll'l and is not very different from recent data from FERMI l lFermi Collaborationll2009bt) : 
the second and the third row are the azimuthal angle and polar angle of radiation source in magnetic polar coordinate, respectively; the 
fourth to the eighth row are for the r^j^ = 1 radiation altitude lower bounds for e = 5 GeV, e = 10 GeV, e = 25 Ge V, e = 100 GeV, 
and £ = 1 TeV, re spectively. During the calculation, we t ake Bq = 3.7 x 10 Gauss jManchester et aLllioOsI) . C, = 60° JWeisskopf et all 
l200Cl : IWan3l2003f) , and p = 0.033 (.Manchcster et al.„2do3) s as Crab pulsar parameters. 



In particular, we calculate the altitudes for last trapping surface, at which the photon generated is optical thick (r = 1) 
due to 7-B processes both analytically 12.11) and numerically ( H2.2p . We regard such altitudes as the lower bounds for 7-ray 
radiation altitude. The analytical estimation and numerical results get good agreement. The numerical calculation suggest 
slightly higher value (20%) for rtr than analytical calculation, due to approximation for analytical method. 

The altitude lower bounds are calculated for different longitudes. For Crab pulsar, the altitude lower bound is about 100 
km, 150 km,250 km, 400 km, and 800 km, for 5 GeV, 10 GeV, 25 GeV, 100 GeV, and 1 TeV7-ray radiation respectively for 
various inclination angle. Note that the rtr is larger around A(/> = 0. This is a geometrical effect. Because the curvature radius 
of magnetic field lines are smaller around Acj> = 0, the transverse magnetic field intensity _Bx and the absorption coefficient 
K, become smaller. Although beyond the measured cut-off energy, 10 GeV pulsed emission is confirmed, 150 k iu radiation 
altitude lower bounds could falsify the classical polar cap model for puslar's 7-ray radiation (IDauehe rtv fc H ardina, 1994 1. 
However it can not discriminat e present pulsar 7-ray radiation model inclu ding slop gap model ( Hardjng ^^aklJg OOSl) . oute r 



gap model (jCheng et al.lll986: 



late present pulsar 7-ray radiation model mclu dmg slop gap model (|JrLarding et aLllzUUell. oute i 
s3Q; Ichiang fc Romani 1994 ; Cheng et al. 2000 h and inner annular gap model ( Qiao et aD 2004ll 



We c an neither exclude possibility for a high altitude model producing a low energy cut-off as indicated in (jHarding et al 
20081 ) ■ since such cut-off may be due to intrinsic physical conditions. 



We see that the lower bounds rtr is not sensitive to Tth, mainly due to the logarithmic dependence of rth, as we have 
already seen in Eq. [l] Thus the difference between rtr for rth = 1 and rth = 10 is less than 10%. 

The magnetic field components of higher order magnetic momentum decrease much faster than that of the dipole magnetic 
momentum. Therefor for the emission from high altitudes (r/R ~ 10~^), we can thus ignore the high er order magnetic momen 



turn than dipole, although the higher order magnetic momentum are important near the star surface (jRuderman fc Sutherland 
I975I ) 



Our results give different value compared to Baring ( 2004h . The main reason is the difference of t reating photon charac- 
teristic absorption length. A brief comparison is given in the appendix of this paper. It is shown that Baring ( 2004h 's result 
under-estimates the characteristic absorption length, such that it needs stronger magnetic field to achieve absorption, which 
under estimate the rtr- 
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The cut-off is not a good radiation location indicator. Fir stly, it is not only sensitive to data quality but also sensitive 

Secondly, there are still pulsed emission 



to the data reduction processes, i.e. how the cut-off is fitted (jLopez et al. 



2OO9I ). 



photons coming from the cut-off tails. Thirdly, one needs to model the initial photon population to make reasonable result 
for cut-off energy. Thus we only use the observation fact that photons of energy e are detected, then we discuss where is 
the lowest possible position (i.e. rtr), above which the photons can just propagate freely without strong absorption by pulsar 
magnetic field. Such rtr is well defined. Future measurement of cut-off energy will not change the conclusion here, as far as 
10 GeV pulse emission of the Crab pulsar has been confirmed, the radiation altitude of 150 km we calculated for the lower 

bounds holds. 

Seven pul sars were detected in 7-ray bands previously ([Thompson 2008), and more than 40 7-ray pulsars are just found 



recently e.g. jFermi CollaborationI boOQal : IPeTlizzoni et al]|2009l 'l. FERMI ( Harding fc Fermi LAT Collaboration! |2009| '1 and 



AGILE are now seeking for more 7-ray pulsars. If young 7-ray pulsars with stronger magnetic field are observed at higher 
energy (e.g. 1 Tev) in the future, there would be a great challenge to the slop gap model and annular gap model. Future 7-ray 
pulsar searching and the follow up observations using ground based Cherenkov telescope ares expected to test present 7-ray 
radiation models for pulsars. 

We are grateful to C.K. Chou for reading the paper. We are also grateful to L. Guillemot and A. Jessner for reading the 
paper and giving valuable comments improving the paper. This work was supported by NSFC (10833003, 10778611, 10821061) 
and the National Basic Research Program of China (Grant 2009CB824800). 
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APPENDIX A: COMPARISON BETWEEN BARING'S RESULTS AND RESULTS IN THIS PAPER 



This paper answers a different questi on compared w ith Barin j 1 2004 ). We are discussing the altitude lower bound for observed 
photons with specific energy, while Barinj ( 2004) di s cusse d the spectra cut-off for collective pho tons. Neverth eless it is 
interesting to see the difference of the results in iBarind |2004l Vs and here. It turns out that we will get lBarin3 (j2004l ')'s results, 
if the characteristic absorption length Ac is taken to be constant, which corresponds to the case where radiation region is near 
the star surface. 

Following Eq. 19 of R uderman and S utherland (1975), the y — B absorption criteria is eMavB±/Bci > Xi where Such 
criteria is also adopted bv iBaring (|2004l ). Due to the curvature of magnetic field, we have B± ~ B sin ~ BXc/rcur ~ 
Bo^^^- So the absorption criteria is 



EMcV 



Bo_ (Roy 2^ ^ 



X 



Substitute rcur = | 



. 27r 



1/2 



which is just 



£Max,[GoV] 



B 



10 



'o \RoJ 



we have 
3 (pr,c/27r)i/2 



Ac 



, ol/2 



/ \ O.lBcr f r,\ 



If we take Ac/x ~ 17 

1/2 



; just get iBarine I 2004 )'s results 
Be 



1,0.1 



GeV. 



(Al) 



(A2) 



(A3) 



(A4) 



If we take Ac = Vs, we go back to Eq. U The reason for making such choice is illustrated in Fig. lAll There are two 
comparative factors dominant the tendency of absorption coefficient k. Firstly, the impact angle 9i between photon direction 
and local magnetic field grows with A then saturated to a limited value due to pure geometrical relation. Secondly, magnetic 
field intensity decrease approximately following (rs + A)~"^ due to dipole field configuration. Thus the perpendicular magnetic 
field strength B± ~ B6i follows relation A/(rs + A)~'^. The absorption coefficient k depends on B± exponentially. So 6i,B and 
K follow the dashed, dot-dashed and solid curve in Fig. I All respectivelv. Clearly the effective absorption length takes value of r. 
Precisely speaking, there must be some numerical fact 7] (due to integration and geometry), where — 1, such that Ac — ryrg. 
In practical, this r] can be regarded as been accounted and put into the threshold optical depth rth. Due to the dependence 
of Ttr on Tth is logarithmic (see Eq.|4]), such correction of rj is not important. This is why the analytical estimation is concord 
with the numerical calculation developed in section [221 

Clearly at high altitude, the characteristic absorption length Ac is no long a constant and must be treated as a f unction 
of radiation altitude. Because of the implicit assumption of a constant characteristic absorption length, BarintJ 1 2004 ) results 
under-estimate the characteristic absorption length Ac so does the optical depth. This leads to smaller vti as we have seen for 
high altitude radiation. We have used Ac ~ as argued in this paper to derive Eq.Hl which is lat er checked by our numerical 
results. For low energy photons (a few GeV), our results agree with Baring fc Hardiii3 ( 2001 ) because the characteristic 
absorption length is about the size of star. 



This paper has been typeset from a T^j]X/ ffl^^X file prepared by the author. 
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Figure Al. The illustration for the relation between 7-B impact angle 6[, magnetic field intensity B, absorption coefliicient re and photon 
propagating length A. 6i,B and re follow the dashed, dot-dashed and solid curve respectively when A grows from to large values. 



